Evolution of Co/Cu multilayer conductivity during growth: 

ab intio study 
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We present ab-initio calculations for the in plane conductivity of Co/Cu multilayer slabs. The 
electronic structure of the multilayer slabs is calculated by means of density functional theory within 
a screened KKR scheme. Transport properties are described using the Boltzmann equation in 
relaxation time approximation. We study the change of the conductivity during growth of the 
multilayer, and we can reproduce the anomalous, non Ohmic, behavior observed experimentally in 
several multilayer systems. Our results show that this behavior can be explained in terms of the 
electronic structure of the slab only. No extra assumption for the scattering at the interfaces is 
necessary. The connection of electronic structure and conductivity during layer-by-layer growth is 
elucidated by analyzing the layer-projected conductivities. 
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I. INTRODUCTION 



The understanding of the size-dependent conductiv- 
ity in thin-film multilayers is still an open question. The 
conductivity of various magnetic/nonmagnetic multilayer 
systemsEfTj has been measured as a function of the thick- 
ness of the deposited layer. These electrical measure- 
ments were made simultaneously by in-situ conductance 
monitoring, and the experimental data show character- 
istic features for all investigated multilayers. First of 
all, the conductance increases with increasing multilayer 
thickness. Second, the slope of the conductance increase 
is different for the different metal layers, due to the corre- 
sponding residual resistivity of the metals. Third, there 
is a characteristic drop in the total conductance as soon 
as the magnetic layer is added on top of the nonmagnetic 
layer. This behavior is quite general and was originally 
related to additional scattering in the interface region due 
to intermixing and disorder. 

Several models have been developed to understand the 
origin of this behavior. Most of them have assumed 
free-electron behavior within the constitucplpLayersLrEl, 
extending the Fuchs-Sondhcimcr approachEjO. Disor- 
der at the interface was introduced. A quantum me- 
chanical description including sipfaee scattering was pre- 
sented bji-.Fishman and CaleckilI3lI3 and by Trivedi and 
AshcrofO. A more realistic description of the multilayer 
was based on a—tight binding model including disorder 
at the interfaced. Ab initio electronic structure calcu- 
lations of magxLctic. multilayers have been performed by 
several groupsLjIij. It was shown that the electronic 
structure of the multilayer plays a very important role 
for the understanding a£..khe transport properties and 
giant magnetoresistanceE9EII. None of the calculations 
were performed for the in-situ conductance-monitoring 
considering the incremental conductance contributed by 
each atomic layer. The aim of this work is the under- 



standing of the size dependence and the evolution of the 
Co/Cu multilayer conductivity during the growth, start- 
ing from ab-initio electronic structure calculations. 



II. COMPUTATIONAL APPROACH 



The calculations are based on Density Functional 
Theory (DFT) in the Local Density Approxima- 
tion (LDA) and a recently implemented version of 
a Screened-, Kprringa-Kohn-Rostoker Green's Function 
mcthodli§l23i2a. For layered systems the advantage of 
this transformation is the N-scaling behavior (with N the 
number of layers in the system) which enables us to treat 
systems with a large number of atoms in the unit cell 
and arbitrary two-dimensional or three-dimensional cell 
periodicity. The potentials were assumed to be spheri- 
cally symmetric inside the Wigner-Seitz sphere as in the 
Atomic Sphere Approximation (ASA) but a multipole 
expansion of the charge density has been taken into ac- 
count up to Imax = 6 (angular momenta up to l rnax = 3 
have been used for the wave functions). For the ex- 
change correlation functionals the expressions given by 
Vosko, Wilk and NusairEj were used. Brillouin zone (BZ) 
integrations hajie been performed by means of special 
points methodsE3. To simulate the in situ conductance- 
monitoring of the Co/Cu multilayer during growth we 
model the system by the following assumptions. We con- 
sider a monolayer growth along the (001) direction. Co 
and Cu are assumed to occupy an ideal fee lattice with 
a lattice constant of 6.76 a.u., without any lattice mis- 
match or distortion at the interfaces and surfaces, while 
no intermixing at the interfaces is allowed. 

In the experimental setup an insulator or semiconduc- 
tor buffer or target material is used, this is simulated by 
vacuum in our calculation. That means a free standing 
slab with 2D periodicity in plane is considered. For the 
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calculations we started from a Cu slab of 2 monolayers 
(ML) thickness. Cu was added layer by layer until a 
thickness of 5 ML which corresponds to about 9A. We 
continue by adding Co up to the same thickness. Thus 
the cycle was repeated. The short-hand notation for the 
considered systems is Cu ni /Co n2 /Cu n3 ... denoting the 
individual layer thicknesses in monolayers starting from 
the bottom Cu layer . Every multilayer corresponds to 
an experimental film after deposition of a complete ide- 
ally flat atomic layer. The electronic structure was de- 
termined self-consistently for each system under consid- 
eration. Concerning the magnetic configuration of the 
mutlilaycr slabs we assumed parallel alignment of the Co 
layer magnetization, which is energetically preferred in 
Co/Cu multilayersEa with 5 ML Cu thickness. 



d 3 i 



l*£(Ri 



(3) 



ASA 



and ^(Ri + r) are normalized wavefunctions of the slab. 
Ri are the basis vectors, denoting the different atomic 
layers, while the real space integration in Eq. (|) is re- 
stricted to the atomic sphere. Thus the total conductivity 
is just the sum over all the layers N, of the slab and spin 
directions 
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IV. ULTRATHIN CU FILM 



III. TRANSPORT THEORY 

To calculate the electrical conductivity of the consid- 
ered layered structures we use the semiclassic al i B pltz- 
mann theory in relaxation time approximation! 16 }! 29 !. A 
spin-independent isotropic relaxation time r is assumed. 
This is of course a rough approximation which concen- 
trates on the anisotropic electronic structure of the sys- 
tem but, as our results show, it is sufficient for the 
explanation of the effect. Moreover, considering the 
anisotropic scattering would change our results only 
quantitatively, and not qualitatively and thus it serves 
our current purpose of understanding the conductivity 
of thin multilayers based on the electronic structure of 
ideal systems. Neglecting spin-flip scattering Mott's two 
current model can be uscdE2l. With these simplifications 
the in-plane conductivity a is given by a Fermi surface in- 
tegral over_the in-plane component of the Fermi velocity 
i£ squaredEil 

ff =ii T E j d 2 \,5{El-E F )vi 2 . (1) 

£7£ are the energy eigenvalues and E F is the Fermi level, 
s denotes the spin directions. The k-integration is per- 
formed over the two-dimensional Brillouin zone of the 
slab. The Fermi surface reduces to a set of closed lines 
in the plane perpendicular to the growth direction . The 
result of Eq. ([!]) is a slab conductivity per unit cell for a 
two-dimensional system. If wc consider only one mono- 
layer then we obtain a conductivity per atom which com- 
pares directly to the bulk conductivity. 

To analyse the spatial distribution of the current a lo- 
cal conductivity is introduced by projecting the conduc- 
tivity on each atomic layer i 

*? = J d 2 k 6(EZ - E F )vi 2 P^{i) (2) 

where 



Following the ideas of Fuchs and SondhcimcrEJii 1 ] the 
conductivity of a simple metallic film goes to zero with 
decreasing film thickness D, since in the limit of thin 
films: a ~ reln(l/K), where k = D/l, and I denotes the 
mean free path. Considering no diffusive surface scatter- 
ing we obtain the bulk conductivity. In a realistic quan- 
tum mechanical calculation where full band structure ef- 
fects are present the situation is more complicated. 
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FIG. 1. Film conductivity (dots) and conductivity 
per atom (crosses) as a function of Cu layer thickness in 
arbitrary units. 

In Fig. [l] wc present the calculated conductivity of ul- 
trathin Cu slabs with thicknesses varying from two to five 
monolayers. The total slab conductivity increases almost 
linearly with the number of Cu monolayers. The conduc- 
tivity per atom saturates to a constant value after 4-5 
ML. It is obvious that our calculated conductivity does 
not extrapolate to zero as predicted by the Sondhcimcr 
model. This is a direct consequence of the quantum me- 
chanical treatment of the surface. In the Sondhcimcr 
model the surface is described by the classical electron 
distribution function. In the present calculation the full 
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quantum mechanical Hamiltonian for the surface is used. 
As a result new types of eigenstates like surface states 
could appear which are not included in the Sondhcimer 
model. Due to the formation of surface states with a 
high in-plane group velocity the conductivity per atom 
(Fig. [I]) is increased for very thin Cu slabs. It has to be 
mentioned that the constant relaxation time that is used 
for all states at the Fermi level might be inaccurate to 
obtain quantitative information. 



V. CONDUCTIVITY OF CU/CO SLABS 
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FIG. 2. Evolution of total and spin projected DOS 
[states/Ryd] at E F and total and spin projected conduc- 
tivity [arbitrary units] during Co/Cu multilayer growth. 



time was assumed, and the difference in the spin chan- 
nels comes solely from the contrast of electronic structure 
in the two spin channels while differences in scattering 
which would require an anisotropic r are ignored in this 
approximation. 

It is interesting to compare the evolution of the conduc- 
tivity with the changes of the slab DOS at the Fermi level, 
shown in Fig. |^. For the majority channel we observe a 
linear increase with thickness since the local DOS at E F 
at the Co and Cu sites is very similar for the majority 
channel and corresponds to the free electron like states 
above the filled 3d bands. The behavior of the minority 
DOS reflects the difference in the electronic structure of 
the two materials. The DOS is much larger at the Co 
sites in comparison to Cu, this causes the strong increase 
of the total DOS during the Co growth. The local DOS 
at the Co surface or at a Co/Cu interface is enhanced in 
comparison to the bulk Co, and the total DOS is strongly 
increased by adding the first two Co layers. Moreover, 
quantum confinement effects are fully included in our cal- 
culation and cause the reduced increase when depositing 
the third Co layer, and the small drop when the single 
Cu layer terminates the slab. 

Let us now concentrate on the drop of the conductiv- 
ity when adding one monolayer Co on the Cu film. This 
behavior is accompanied by a strong increase in the total 
DOS since the Fermi level crosses the Co minority band. 
The spin projected conductivities demonstrate that the 
drop occurs in the minority channel. By means of an 
averaged spin-dependent Fermi velocity v B 



d 2 kS(E^- E F )vl 2 /n s {E F ) 



(5) 



one can rewrite Eq. (Q) into 



Our calculated conductivity of the Cu/Co slabs is 
shown in Fig. |^, together with the DOS at the Fermi 
level. With our calculation we are able to reproduce the 
experimentally observed behavior. Our results show an 
increase with increasing slab thickness, while there is a 
drop in the slab conductivity of about 20% when the 
first Co monolayer is deposited on Cu. By adding more 
Co layers a increases further, and we observe no signif- 
icant effect when we start depositing Cu again. A drop 
is again seen when we repeat the Co deposition on top 
of Cu. The spin decomposition shows that for majority 
electrons we have a monotonous increase in er, while the 
peculiar, non-Ohmic conductivity variation stems from 
the minority electrons. Here we point out that no specific 
surface scattering mechanism is included in the calcula- 
tion. Moreover, a constant spin-independent relaxation 



n s (E F ) = J d 2 k 6(E£ - E F ) is the total DOS at Fermi 
level. The expression in Eq. (^|) is useful to understand 
the microscopic origin of the conductivity drop. Follow- 
ing Eq. ^ the reduction of the conductivity is caused by 
a strong reduction of the averaged Fermi velocity in the 
minority channel since the DOS is continously increas- 
ing. For the majority band we have a linear increase and 
no difference is seen when going from Cu to Co and vice 
versa since the Co d states are below the Fermi level. 
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FIG. 3. Evolution of the local conductivity of the 
minority channel (Eq. ^) [arbitrary units] during film 
growth. 




To elucidate this Fermi velocity reduction the local 
conductivity of the minority channel (Eq. ||) is shown in 
Fig. ||, for nearly all slabs under consideration. For the 
pure Cu slabs all Cu layers contribute equally to the con- 
ductivity, and we observe only a small reduction at the 
surface layers. The situation changes drastically when 
we consider a single monolayer of Co on the Cu slab. 
As we can sec from Fig. || (third panel from bottom), 
the major contribution in the conductivity comes from 
the Co layers, the Cu contributions however, are signifi- 
cantly decreased by almost a factor of four in comparison 
to the Cu slabs without Co coverage. According to Eq. 
(||) , this corresponds to a reduction of the effective Fermi 
velocity by a factor of two and the total slab conductivity 
is decreased (Fig. ||). 

This reduction of the Fermi velocity stems from a dras- 
tic change of the Fermi surface topology when adding a 
single Co layer on top of Cu as shown in Fig. |[ The 
deposition of a simple Co layer destroys the free electron 
like Cu-slab Fermi surface in the minority band while 
the majority Fermi surface is kept intact. As a result we 
end up with a complicated minority Fermi surface dom- 
inated by Co d states of low velocity (Fig. ^). Despite 
the low velocity of the Co d-states, Co layers have a high 
a because of the high DOS. Moreover, Co surface and 
interface states produce a further enhancement. 




Cu§ Cu§Co\ 
FIG. 4. Fermi surface of Cms and Cu$Coi slabs 
for minority (upper panel) and majority (lower panel) 
bands. Note that for 2d slabs of finite thickness the Fermi 
surface reduces to a set of lines. 



The constant reduction of the Cu contribution implies 
also the experimentally found proportionality of the coft- 
ductivity drop in dependence on the Cu layer thicknesscl. 
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The conducting states have a large probability ampli- 
tude at the surface, which would give an amplification 
of the conductivity drop due to the enhanced scattering 
cross sectiojpEj. This is in line with numerical results 
of TsymbaLl assuming a larger potential disorder at the 
Co/Cu interface in comparison to the bulk materials. 

The addition of any extra Cu layer on top of the Co 
layer is accompanied by a similar conductivity contribu- 
tion like in the first Cu stack. In the minority chan- 
nel there is only a marginal decrease when we add one 
monolayer Cu on top of Co, this is compensated by the 
linear increase in the majority channel, so that almost 
no effect is observed in the total a. With increasing Cu 
layer thickness a symmetric profile is formed (Fig. sec- 
ond panel from top). The total conductivity increases 
monotonously. The first Co layer of the second stack 
causes again a decrease of the underlying Cu layer con- 
tributions and a drop of the total a. The combination 
of all details in the local conductivity discussed above 
leads to the observed behavior as a function of the film 
thickness (Fig. |J), decreasing conductivity adding Co on 
Cu, and a linear increase when adding Cu on Co. The 
Cu thickness we consider is much smaller compared with 
the recent experiment of Bailey et aid. but as we can see 
from Figs, [I] and |[ a constant conductivity per atom is 
established after 4-5 ML and larger thickness would only 
cause a monotonic increase in a. As mentioned above 
our calculations are performed for a ferromagnetic con- 
figuration of the Co magnetization which is favored by 
the interlayer distance of 5 ML of Cu. The Cu thick- 
ness, of 20A (11 ML) used in the experiment by Bailey et 
al.E3 and Urbaniak et al.ul corresponds to the second anti- 
ferromagnetic maximum of interlayer exchange coupling 
in Co/Cu and Py/Cu multilayers, respectively. However 
the actual magnetic configuration of the multilayer slab 
under consideration was not investigated in the experi- 
ments. 

The microscopic picture discussed above could be ap- 
plied to other systems with similar electronic structure 
like Py/CuO, Fc/Agta, and Cr/AuH. Depending on the 
electronic structure of the particular system the role of 
minority and majority electrons can be changed. The 
only system that has to be checked in more detail is 
Cu/Agp. Cu and Ag are isovalent and the bandstruc- 
tures arc quite similar. The conductivity drop in this 
system might be related to alloying only. 



VI. CONCLUSIONS 

Based on ab-initio calculations of the conductivity of 
Cu/Co slabs we have shown that the variation during 
growth can be explained by the changes of the electronic 
structure as a function of the film thickness. In general 
the total conductivity of a multilayer slab increases with 
increasing thickness. A pronounced conductivity drop is 
obtained when the first Co layer is deposited on top of the 



Cu. Due to interaction of the Co d-electrons with the Cu 
s-electrons the character of the minority Fermi surface is 
changed from sp to d-like. As a consequence the local 
conductivity contributions of the Cu layers are reduced. 
Since the total DOS at the Fermi level is increasing dur- 
ing film growth the effect is related to reduced Fermi 
velocities. The microscopic picture can be generalized to 
any multilayer consisting of noble and transition metal 
layers. This effect is not in contradiction to the con- 
duction dron-caused by interface scattering discussed by 
Bailey et al£3 but can be combined with it. Both ef- 
fects exist and interface scattering can even amplify the 
bandstructure effects. 



1 W.E. Bailey, S.X. Wang, and E.Y. Tsymbal, Phys. Rev. B 
61, 1330 (2000). 

2 H. Briickl, PhD Thesis, Universitat Regensburg (1992), 

3 H. Eckl, G. Reiss, H. Briickl, and H. Hoffmann, J. Appl. 
Phys. 75, 362 (1994). 

4 S. Fahler, M. Weisheit, and H.-U. Krebs, Mat. Res. Soc. 
Proc. 502, 139 (1998). 

5 M. Urbaniak, T. Lucihski, and F. Stobiecki, J. Mag. Mag. 
Mat. 190, 187 (1998). 

e S. Fahler, M. Weisheit, K. Sturm, and H.-U. Krebs, Appl. 
Surf. Sci. 154-155, 419 (2000). 

7 G. Palasantzas and J. Barnas, Phys. Rev. B 56, 7726 
(1997). 

8 W.A. Harrison, Phys. Rev. B 61, 7766 (2000). 

9 G. Palasantzas, Y.-P. Zhao, G.-C. Wang, T.-M. Lu, J. 
Barnas, and J. Th. M. De Hosson, Phys. Rev. B 61, 11109 
(2000). 

10 K. Fuchs, Proc. Royal Soc. A 34, 100 (1938). 

11 E. H. Sondheimer, Adv. Phys. 1, 1 (1952). 

12 G. Fishman and D. Calecki, Phys. Rev. Lett. 62, 1302 
(1989). 

13 G. Fishman and D. Calecki, Phys. Rev. B 43, 11581 (1991). 

14 N. Trivedi and N. W. Ashcroft, Phys. Rev. B 38, 12298 
(1988). 

15 K.M. Schep, P.J. Kelly, and G.E.W. Bauer, Phys. Rev. 

Lett. 74, 586 (1995). 
16 P. Zahn, I. Mertig, M. Richter, H. Eschrig, Phys. Rev. Lett. 

75, 2996 (1995). 
17 W.H. Butler, X.-G. Zhang, D.M.C. Nicholson, T.C. 

Schulthess, and J.M. MacLaren, Phys. Rev. Lett. 76, 3216 

(1996). 

18 P. Zahn, J. Binder, I. Mertig, R. Zeller, and P.H. Dederichs, 
Phys. Rev. Lett. 80, 4309 (1998). 

19 C. Blaas, P. Weinberger, L. Szunyogh, P.M. Levy, and C. 
Sommers, Phys. Rev. B 60, 492 (1999). 

20 G. Binasch, P. Griinberg, F. Saurenbach, and W. Zinn, 
Phys. Rev. B 39, 4828 (1989). 

21 M.N. Baibich, J.M. Broto, A. Fert, F. Nguyen Van Dau, 
F. Petroff, P. Etienne, G. Creuzet, A. Friederich, and J. 
Chazelas, Phys. Rev. Lett. 61, 2472 (1988) 



5 



O.K. Anderson, A.V. Postnikov, and S.Yu. Savrasov, Mat. 
Res. Soc. Proc. 253, 37 (1992). 

L. Szunyogh, B. Ujfalussy, P. Weinberger, and J. Kollar, 
Phys. Rev. B 49, 2721 (1994). 

R. Zeller, P.H. Dederichs, B. Ujfalussy, L. Szunyogh, and 
P. Weinberger, Phys. Rev. B 52, 8807 (1995). 
K. Wildberger, R. Zeller, and P.H. Dederichs, Phys. Rev. 
B 55, 10074 (1997). 

S. H. Vosko, L. Wilk, and M. Nusair, Can. J. Phys. 58, 
1200 (1980). 



H. J. Monkhorst and J.D. Pack, Phys. Rev. B 13, 5188 
(1976). 

P. Lang, L. Nordstrom, K. Wildberger, R. Zeller, P.H. Ded- 
erichs, and T. Hoshino, Phys. Rev. B 53, 9092 (1996). 

I. Mertig, Rep. Prog. Phys. 62, 237 (1999). 
N.C. Mott, Adv. Phys. 13, 325 (1964) 

P. Zahn, I. Mertig, R. Zeller, and P.H. Dederichs, Mat. Res. 
Soc. Proc. 475, 525 (1997). 

W.E. Bailey, C. Fery, K. Yamada, and S.X. Wang, J. Appl. 
Phys. 85, 7345 (1999). 



G 



